Neutron scattering as a probe of the Fe-pnicitide superconducting gap 
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Inelastic neutron scattering provides a probe for studying the spin and momentum structure of 
the superconducting gap. Here, using a two-orbital model for the Fe-pnicitide superconductors and 
an RPA-BCS approximation for the dynamic spin susceptibility, we explore the scattering response 
for various gaps that have been proposed. 



Introduction - Recent neutron scattering experiments 
have shown that LaOFeAs undergoes a structural distor- 
tion below ~150K, which is then followed at ~137K by 
the onset of long-range spin density wave (SDW) order 
with a wave vector q = (0.5, 0.5, 0.5)7r/a When it 

is doped with F, both the structural distortion and the 
magnetic order are suppressed and the system becomes 
superconducting []]]. When La is replaced by Sm, su- 
perconducting transition temperatures of 55K have been 
reported [IJ. Thus it is natural to believe that the 
Fe-pnicitide superconductors have an electronic pairing 
mechanism and a variety of unconventional ga p st ruc- 
tures have already been proposed 

Here we explore how inelastic neutron scattering in the 
superconducting state can provide information on which 
of the gap structures actually occurs. 

Bandstructure calculations for doped LaOFeAs give a 
Fermi surface for the 2-Fe/cell Brillouin zone which con- 
sists of two hole cylinders around the r_-point and two- 
electron cylinders around the M-point [J, la, LL2J • When 
this is folded out into the larger Brillouin zone which 
is associated with a square lattice of Fe sites having 1- 
Fe/cell, one gets a similar Fermi surface to that shown 
in Fig. la [4j. In this case the SDW would be associ- 
ated with q = (w, 0) or (0, 7r). Here we will use a simple 
2-orbital per site tight-binding model, which has been 
parametrized to give the Fermi surface shown in Fig. 1 
hj. We then assume that the spin susceptibility in the 
superconducting state can be modeled by an RPA-BCS 
form and proceed to explore the structure of the inelastic 
scattering in the superconducting state for two q values 
and various gaps. As one knows from the cuprate prob- 
lem, the occurance of resonances in the neutron scatter- 
ing depend through the BCS coherence factors on the 
relative signs of the gap on different parts of the Fermi 
surface which are separated by q. For the present case, 
in which there are multiple Fermi surfaces, there is a va- 
riety of ways in which resonances can occur and provide 
information on the gap structure. 

In the following, we first give a brief review of the 
model and then outline the RPA-BCS calculation of the 
spin susceptibility. This is similar to various approxi- 
mations used for both the cuprates [14J and Sr2Ru04 



|15j |. We calculate the normal RPA spin susceptibility 
and then examine the RPA-BCS response in the super- 
conducting state for both singlet and triplet gap func- 
tions. A related study based on a four-band model was 
reported in Ref. [l6| for the case of a singlet gap. How- 
ever, these authors neglected to take into account the 
matrix elements which relate the band operators to the 
orbital operators. The singlet gap functions we will use 
correspond to low order lattice harmonic representations 
of the sign-reversed s-wave gap proposed by Mazin et al. 



A s *(k) = Ao cos k x cos k y 
and an extended s-wave gap 

A xs (k) = Ao/2(cosfc :E +cosfc y ) 



(1) 



(2) 



We will aso determine the spin-flip Imx^ (q, oS) and non- 
spin- flip Imx Z z(q, w) response for various p-wave triplets. 
In this case, 



A a/3 (k) = d(k) • aia 2 



aft 



(3) 



with d perpendicular to the Fe plane and d z (k) — A(k). 
Here we will explore A(k) = sinfe^, sin 2k x and sinfc x + 
i sin k y . 

We conclude with a summary of what one can expect 
to learn about the gap symmetry from inelastic neutron 
scattering in the superconducting state. 

Model - Bandstructure calculations for doped 
LaOFeAs show that the low energy states near the Fermi 
energy have dominant 3d Fe character and various multi- 
orbital tight-binding fits have been proposed. Here we 
will use a minimal, 2-orbital "d xz -d yz v per site tight- 
binding model with parameters chosen to give the Fermi 
surfaces shown in Fig. la. This model has the virtue of 
simplicity while qualitatively capturing the shapes of the 
bandstructure Fermi surfaces and the relationship be- 
tween the band operators to the orbital operators. This 
latter feature is important since it is the variation of the 
gaps on the Fermi surfaces that determine, through the 
BCS coherence factors, the inelastic neutron scattering 
response. While the magnitude of the response depends 
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upon the Fermi velocities on the Fermi surfaces, which 
are not well reproduced by the 2-orbital model, the occu- 
rance or non-occurance of resonant features is determined 
by the k-dependence of the gap and the Femri surfaces. 

As described in Ref. [131] . our minimal model consists 
of a square two-dimensional lattice with degenerate " d xz " 
and"(iy Z " orbitals on each site. One-electron hopping pa- 
rameters U are introduced which provide near-neighbor 
<j(t\) and 7r(£ 2 ) couplings between similar orbitals, as well 
as a next near-neighbor coupling £3. In addition, there 
is a next-near-neighbor coupling t$ which hybridizes d xz 
with d yz . The resulting tight-binding Hamiltonian can 
be written as 



Ho 



with 



ko- 



^(k) [( e+ (k) - fi)I + £_(k)r 3 + ex»(k)n] V CT 



d ytT (k) 



(4) 
(5) 



Here Tj are the usual Pauli matrices and 
i[ eiB (k)±e v (k)], 



2^2 cos fc„ — 4ts cos fc T cos fc„ 



e±(k) 

£a;(k) = —2tl COS fcj; — Zt2 COS ftj, — It 3 Oud n, x <^wo . 

E y (k) — — 2t2 cos k x — 2t\ cos k y — 4i 3 cos fca, cos k y , 

e xl/ (k) = — 4^4 sin fc s sinfcy . (6) 



Introducing the band-operator j u(T (k), such that for r 
x or y 



Vw(k) 



a^(k)7 1/er (k) 



(7) 



with 



ai(k) - o»(k) = sgn(e xy (k)) 



\J 2 2^(k)+4„(k) 



(k) 



<(k) 



-«-(k) 



\ 



e_(k) 



2^(k)+ £ yk) 
the tight-binding Hamiltonian becomes 

ff = E ^(k)7L(k)7„(k). 

Here the band-energies are 



E ± (k) = £+ (k) ± J el (k) + £ 2 ( k ) - M . 



(8) 



(9) 



(10) 



We choose the hopping parameters t\ = —1, t 2 = 1.3, 
<3 = <4 = —0.85 and chemical potential /i = 1.45. The ol\ 
and «2 Fermi surfaces in Fig. la correspond to E- (kp) = 
0, while the /?i and /?2 surfaces correspond to E + (kp) = 
0. 



Inelastic neutron scattering - The physical spin sus- 
ceptibility 



' (q, iu) m ) = E Xrs(<l, «^ m ) 



(11) 



is calculated from the orbital dependent spin susceptibil- 
ity defined as 

^M,iw m ) = J ' dT(T T Sl(q,r)S s 3 (-q,0)). (12) 



Here, r, s — x,y label the orbital indices, and 5f(q) = 
i X)k' ! / ; ra(k + q)er^Vv/3(k) is the ith-component of the 
spin operator for orbital r. In the BCS framework, one 
obtains for the BCS orbital dependent spin susceptibility 

xt.(<h %) = E M " ' ( k ' q) x ( 13 ) 

k,n 

{GL(k + q)G£>) + F#(-k - q)F^(k)} . 

Here we used fc = (k, u> n ) and g = (q, ui m ) an d f, 
■ . — are the eigenvalues of the bands. The normal and 
anomalous Green's functions for band v are given by 



iu m + E u (k) 



A a/J (k) 



(14) 



with £„(k) = y/EZ(k) + |A(k)| 2 . The hybridization be- 
tween the bands is reflected in the matrix elements 

M r 7'(k,q) =<(k + q)a^(k)<,(k)<(k + q). (15) 

One can then obtain the BCS spin susceptibilty 

Xo (q, oj) — J2 r ,s Xo^rsCl) u ) on tne rea l frequency axis 
from an analytical continuation of Eq. (fT5|) . We then 
use the RPA to take into account the effect of the on- 
site intra-orbital Coulomb interactions U [l]| ■ The RPA 
susceptibility is determined from the matrix equation 



Xrpa(«) = E 



X \q)(l-Tx \q) 



where q — (q, a;), with the interaction vertex 

U 
U 



(16) 



(17) 



Fig. Hb shows the results for the imaginary part of the 
RPA spin susceptiblity Xo"rpa( c 1' w ) m ^ ne normal state 
(A = 0) for q = (tt, 0) and q* = (n/2, 0). Here we used 
U = 3|* x 1, [i = 1.45 and the temperature T = 0.001|ti|. 

For fj, — 1.45, the static RPA spin susceptibility shows 
peaks at q = (n,0) and q* = (7r/2,0). As indicated by 
the dotted and dashed arrows in Fig. la, q = (tt, 0) is a 
nesting vector between a region on the a\ Fermi surface 
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(a) 



P2 \ a 2 ~ 
ai \ Pi 



q=(3t,0)' 



q*=(ic/2,0) 




0.20 



FIG. 1: (a) Fermi surface of the two-orbital model in the large 
1-Fe/cell Brillouin zone. When this is folded down along the 
dashed line, one obtains the Fermi surface of the 2-Fe/cell 
system. The ai/% Fermi surface sheets are hole pockets given 
by _E-(1cf) = and the /3i/2 sheets are electron pockets 
given by E+(Is.f) = 0. The nesting vectors q = (7r,0) and 
q* = (7r/2,0) are indicated by the dotted and dashed lines, 
respectively, (b) RPA spin susceptibility Xrpa^i versus 
frequency in the normal state (A(k) = 0) for q = (71", 0) and 
q*=(7r/2,0). 



tors that enter the spin susceptibility depend upon the 
sign of A(k + q)A(k). For a singlet gap, when this 
is negative, there can be a resonance response at w = 
A(k + q)| + |A(k)|. For a triplet gap, the coherence fac- 
tors are different for Xzz and \^ I n this case, when 

Re(A*(k + q)A(k)) is negative, there can be a reso- 
nance in Xzz but not in \^ Likewise, when Re(A*(k + 

q)A(k)) is positive, X-\ can exhibit a resonance while 

Xzz varies smoothly through w = |A(k + q)| + |A(k)|. 
Thus, the neutron scattering response in the supercon- 
ducting state can provide information on the momentum 
and spin structure of the superconducting gap. Here we 
examine the response for various gaps that have been 
proposed for LaOFeAs. 
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sheet with the Pi Fermi surface sheet, while q* = (tt/2, 0) 
connects several different regions on the Fermi surface. 
The scattering for q* = (w/2,0) is, however, dominated 
by the intra-band excitations on the 0% Fermi surface 
sheet (see Fig. lb), since for the other inter-band pro- 
cesses s xy (k) and e xy (k+q*) are close to zero, and there- 
fore the matrix-elements M"g (k, q) are very small. At 
low frequency Xrpa( < 1' w ) 1S iar ger for q = (tt, 0) than for 
q = (tt/2, 0), since larger regions of the Fermi surface are 
nested for q = (71-, 0). At low frequency the RPA sus- 
ceptibility is strongly enhanced over the unrenormalized 
susceptibility. 

In the superconducting state, the gap A(k) is finite and 
the susceptibility, Eq. (fT3"|) , depends on the symmetry of 
the gap. For singlet pairing one has Fj|(k) = —F^(k) 
and i r fx(— k) = Fjj(k) and hence the in-plane suscep- 
tibility x H = \{x xx + X vv ) 1S equal to the out-of- 
plane susceptibility x zz as m the normal state. For 
the triplet case, however, one has F-fj(k) = F^(k) and 
i*tj,(— k) = — -Fjx(k). Therefore x H an d X zz differ with 
respect to their superconducting coherence factors. 

Results - As is well known, the BCS coherence fac- 



FIG. 2: RPA dynamic spin susceptibility Xr,pa(ci> aj ) versus 
frequency in the spin singlet superconducting state for (a) 
q = (tt, 0), and (b) q* = (tt/2, 0). 

To begin, Fig. 2a shows the imaginary part of the RPA- 
BCS spin susceptibility for various singlet gaps at a mo- 
mentum transfer q = (7r,0). We have modeled a sign- 
reversed s-wave gap A s *(k), i.e. the type proposed by 
Mazin et al. by Eq. fl} and the solid line in Fig. [2k 
shows the expected resonance response associated with 
having A(k + q)A(k) < 0. This behavior can be con- 
trasted with the response found for an extended s-wave 
gap. Here A(k + q)A(k) > for the q = (ir,0) nesting 
vector. 

Fig. 2b shows the results for q* = (7r/2,0). For this 
momentum transfer the scattering is dominated by the 
intra-band process connecting regions on the [3\ Fermi 
surface sheet. For these regions one has k + q* = k, 
and therefore A(k + q*)A(k) > for the singlet gaps. 
Thus, no resonance is found for q* = (7r/2, 0) in this case. 

Similar results for various triplet gaps are shown in 
Fig. [3] for q = (tt, 0). For the p^-wave sinfc^ gap, 
A(k + q)A(k) < and a resonance is seen in x'zzi'i' w ) 
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FIG. 3: RPA spin susceptibility (a) X+-,rpa(cU w ) an d (b) 
Xzz rpa( ( 1i w ) versus frequency in the spin triplet supercon- 
ducting state for q = (ir, 0). 



but not in x+— (<1> W )- One could also consider a sign- 
reversed p-wave modeled by A P 2 X — Aosin2fc a; . Here 
the coherence factor for q = (71", 0) is positive so that the 
resonance appears in x'+-{ ( \t 1jJ )- For the sinfc K + ism.k v 
gap, smk x sva{k x + q x ) < and sin k y sm{k y + q y ) > 
with similar size for the dominant process. Hence, 
Re(A*(k + q)A(k)) is close to zero and one obtains qual- 
itatively similar results for x~\ an d Xzz- 




r 8. 



00 0.05 0.10 0.15 

CO 



0.00 0.05 0.10 0.15 0.20 

CO 



FIG. 4: RPA spin susceptibility (a) X+-,rpa(cL w ) an d (b) 
Xzz rpa( ( 1i w ) versus frequency in the spin triplet supercon- 
ducting state for q* = (tt/2, 0). 



Fig. 4 shows the results for the triplet gaps for q* = 
(7r/2, 0). For the p^-wave smk x gap, the results are sim- 



ilar as in Fig. 3 for q = (tt, 0). The gap changes sign 
under the transformation k — > — k so that the gap has 
opposite signs on the sheets of the /3i Fermi pocket con- 
nected by q* (see Fig. la). Hence, x'zzi^ 1 ^) displays a 
resonance while x"_(q, uj) does not. For the sin2fc x gap, 
the situation is similar and a resonance is found in Xzz, 
but not in X+-- This is opposite to the results in Fig. 3 
for q = (tt, 0). The results for the sinfc x + isink y case 
are almost identical to the results found for the p^-wave 
gap. This is explained by the fact that for the domi- 
nant intra-band scattering on the /3i Fermi surface sheet, 
k y and k y + q y is close to zero, and therefore the sinfcy 
contribution to the p x +ip y gap is an order of magnitude 
smaller than the sink x contribution. 

Conclusion - Using a two-orbital model for the Fe- 
pnicitide superconductors and an RPA-BCS approxima- 
tion for the dynamic spin susceptibiltiy we have explored 
the inelastic scattering response for various gaps that 
have been proposed. As one would expect, we have found 
that the occurance of resonances in the dynamic spin sus- 
ceptibility depends on the relative signs of the gap on the 
parts of the Fermi surface separated by q, and in the case 
of triplet gaps also on whether the in-plane, %"_(q, u>) or 
out-of-plane, x" z (q, a;) components are studied. Specif- 
ically, for the singlet gaps, we have found that x"{ c Lt uj ) 
displays a resonance for a sign-reversed s-wave gap mod- 
eled by Ao cos k x cos k y for q = (tt, 0) which corresponds 
to the antiferromagnetic wave- vector, but no resonance 
for the extended s-wave gap or for q* = (7r/2,0). For 
the triplet case, resonances were found for a p x -wave gap 
in x'lz an d for a sign-reversed p-wave gap in X+- f° r 
q = (tt, 0). For q* = (n/2,0), resonances appeared in 
x'z Z for the p-wave, sign- reversed p-wave andp^+ipy-wave 
gaps. 
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